Birmingham Hip Resurfacing at 20 years.

Aims: The Birmingham Hip Resurfacing (BHR) arthroplasty has been used as a surgical treatment of coxarthrosis since 1997. We present 20-year results of 234 consecutive BHRs performed in our unit. Methods: Between 1999 and 2001, there were 217 patients: 142 males (65.4%), mean age 52 years (18 to 68) who had 234 implants (17 bilateral). They had patient-reported outcome measures collected, imaging (radiograph and ultrasound), and serum metal ion assessment. Survivorship analysis was performed using Kaplan-Meier estimates. Revision for any cause was considered as an endpoint for the analysis. Results: Mean follow-up was 20.9 years (19.3 to 22.4). Registry data revealed that 19 hips (8.1%) had been revised and 26 patients (12%) had died from causes unrelated to the BHR. Among the remaining 189 hips, 61% were available for clinical follow-up at 20 years (n = 115) and 70% of patients had biochemical follow-up (n = 132). The cumulative implant survival rate at 20 years for male patients was 96.5% (95% confidence interval (CI) 93.5 to 99.6), and for female patients 87% (95% CI 79.7 to 94.9). The difference was statistically significant (p = 0.029). The mean Oxford Hip Score, Hip disability and Osteoarthritis Outcome Score, and Forgotten Joint Score were 45 (29 to 48), 89 (43 to 100), and 84 (19 to 100), respectively. The mean scores for each of the five domains of the EuroQol five-dimension three-level questionnaire were 1.2, 1.0, 1.2, 1.3, and 1.1, and mean overall score 82.6 (50 to 100). Ultrasound showed no pseudotumour. Mean cobalt and chromium levels were 32.1 nmol/l (1 to 374) and 45.5 nmol/l (9 to 408), respectively. Conclusion: This study shows that BHRs provide excellent survivorship and functional outcomes in young male patients. At 20 years, soft-tissue imaging and serum metal ion studies suggest that a metal-on-metal resurfacing implant can be well tolerated in a group of young patients.

DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation.

BACKGROUND: Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions. RESULTS: The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression. CONCLUSIONS: In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.

DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation

BACKGROUND: Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions. RESULTS: The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression. CONCLUSIONS: In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.

HDAC6 regulates NF-κB signalling to control chondrocyte IL-1-induced MMP and inflammatory gene expression.

Elevated pro-inflammatory signalling coupled with catabolic metalloproteinase expression is a common feature of arthritis, leading to cartilage damage, deterioration of the joint architecture and the associated pain and immobility. Countering these processes, histone deacetylase inhibitors (HDACi) have been shown to suppress matrix metalloproteinase (MMP) expression, block cytokine-induced signalling and reduce the cartilage degradation in animal models of the arthritis. In order to establish which specific HDACs account for these chondro-protective effects an HDAC1-11 RNAi screen was performed. HDAC6 was required for both the interleukin (IL)-1 induction of MMP expression and pro-inflammatory interleukin expression in chondrocytes, implicating an effect on NF-κB signalling. Depletion of HDAC6 post-transcriptionally up-regulated inhibitor of κB (IκB), prevented the nuclear translocation of NF-κB subunits and down-regulated NF-κB reporter activation. The pharmacological inhibition of HDAC6 reduced MMP expression in chondrocytes and cartilage collagen release. This work highlights the important role of HDAC6 in pro-inflammatory signalling and metalloproteinase gene expression, and identifies a part for HDAC6 in the NF-κB signalling pathway. By confirming the protection of cartilage this work supports the inhibition of HDAC6 as a possible therapeutic strategy in arthritis.

Matrix metalloproteinase-13 is fully activated by neutrophil elastase and inactivates its serpin inhibitor, alpha-1 antitrypsin: Implications for osteoarthritis.

Matrix metalloproteinase-13 (MMP-13) is a uniquely important collagenase that promotes the irreversible destruction of cartilage collagen in osteoarthritis (OA). Collagenase activation is a key control point for cartilage breakdown to occur, yet our understanding of the proteinases involved in this process is limited. Neutrophil elastase (NE) is a well-described proteoglycan-degrading enzyme which is historically associated with inflammatory arthritis, but more recent evidence suggests a potential role in OA. In this study, we investigated the effect of neutrophil elastase on OA cartilage collagen destruction and collagenase activation. Neutrophil elastase induced significant collagen destruction from human OA cartilage ex vivo, in an MMP-dependent manner. In vitro, neutrophil elastase directly and robustly activated pro-MMP-13, and N-terminal sequencing identified cleavage close to the cysteine switch at 72 MKKPR, ultimately resulting in the fully active form with the neo-N terminus of 85 YNVFP. Mole-per-mole, activation was more potent than by MMP-3, a classical collagenase activator. Elastase was detectable in human OA synovial fluid and OA synovia which displayed histologically graded evidence of synovitis. Bioinformatic analyses demonstrated that, compared with other tissues, control cartilage exhibited remarkably high transcript levels of the major elastase inhibitor, (AAT) alpha-1 antitrypsin (gene name SERPINA1), but these were reduced in OA. AAT was located predominantly in superficial cartilage zones, and staining enhanced in regions of cartilage damage. Finally, active MMP-13 specifically inactivated AAT by removal of the serine proteinase cleavage/inhibition site. Taken together, this study identifies elastase as a novel activator of pro-MMP-13 that has relevance for cartilage collagen destruction in OA patients with synovitis.

Osteoporosis, fracture, osteoarthritis & sarcopenia: A systematic review of circulating microRNA association.

Circulating microRNAs (c-miRs) show promise as biomarkers. This systematic review explores their potential association with age-related fracture/osteoporosis (OP), osteoarthritis (OA) and sarcopenia (SP), as well as cross-disease association. Most overlap occurred between OA and OP, suggesting potentially shared microRNA activity. There was little agreement in results across studies. Few reported receiver operating characteristic analysis (ROC) and many identified significant dysregulation in disease, but direction of effect was commonly conflicting. c-miRs with most evidence for consistency in dysregulation included miR-146a, miR-155 and miR-98 for OA (upregulated). Area under the curve (AUC) for miR-146a biomarker performance was AUC 0.92, p = 0.028. miR-125b (AUC 0.76-0.89), miR-100, miR-148a and miR-24 were consistently upregulated in OP. Insufficient evidence exists for c-miRs in SP. Study quality was typically rated intermediate/high risk of bias. Wide study heterogeneity meant meta-analysis was not possible. We provide detailed critique and recommendations for future approaches in c-miR analyses based on this review.

Increased hippocampal excitability in miR-324-null mice.

MicroRNAs are non-coding RNAs that act to downregulate the expression of target genes by translational repression and degradation of messenger RNA molecules. Individual microRNAs have the ability to specifically target a wide array of gene transcripts, therefore allowing each microRNA to play key roles in multiple biological pathways. miR-324 is a microRNA predicted to target thousands of RNA transcripts and is expressed far more highly in the brain than in any other tissue, suggesting that it may play a role in one or multiple neurological pathways. Here we present data from the first global miR-324-null mice, in which increased excitability and interictal discharges were identified in vitro in the hippocampus. RNA sequencing was used to identify differentially expressed genes in miR-324-null mice which may contribute to this increased hippocampal excitability, and 3'UTR luciferase assays and western blotting revealed that two of these, Suox and Cd300lf, are novel direct targets of miR-324. Characterisation of microRNAs that produce an effect on neurological activity, such as miR-324, and identification of the pathways they regulate will allow a better understanding of the processes involved in normal neurological function and in turn may present novel pharmaceutical targets in treating neurological disease.

Dynamic chromatin accessibility landscape changes following interleukin-1 stimulation.

Dynamic modifications of chromatin allow rapid access of the gene regulatory machinery to condensed genomic regions facilitating subsequent gene expression. Inflammatory cytokine stimulation of cells can cause rapid gene expression changes through direct signalling pathway-mediated transcription factor activation and regulatory element binding. Here we used the Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) to assess regions of the genome that are differentially accessible following treatment of cells with interleukin-1 (IL-1). We identified 126,483 open chromatin regions, with 241 regions significantly differentially accessible following stimulation, with 64 and 177 more or less accessible, respectively. These differentially accessible regions predominantly correspond to regions of the genome marked as enhancers. Motif searching identified an overrepresentation of a number of transcription factors, most notably RelA, in the regions becoming more accessible, with analysis of ChIP-seq data confirmed RelA binding to these regions. A significant correlation in differential chromatin accessibility and gene expression was also observed. Functionality in regulating gene expression was confirmed using CRISPR/Cas9 genome-editing to delete regions that became more accessible following stimulation in the genes MMP13, IKBKE and C1QTNF1. These same regions were also accessible for activation using a dCas9-transcriptional activator and showed enhancer activity in a cellular model. Together, these data describe and functionally validate a number of dynamically accessible chromatin regions involved in inflammatory signalling.

OATargets: a knowledge base of genes associated with osteoarthritis joint damage in animals.

OBJECTIVES: To collate the genes experimentally modulated in animal models of osteoarthritis (OA) and compare these data with OA transcriptomics data to identify potential therapeutic targets. METHODS: PubMed searches were conducted to identify publications describing gene modulations in animal models. Analysed gene expression data were retrieved from the SkeletalVis database of analysed skeletal microarray and RNA-Seq expression data. A network diffusion approach was used to predict new genes associated with OA joint damage. RESULTS: A total of 459 genes were identified as having been modulated in animal models of OA, with ageing and post-traumatic (surgical) models the most prominent. Ninety-eight of the 143 genes (69%) genetically modulated more than once had a consistent effect on OA joint damage severity. Several discrepancies between different studies were identified, providing lessons on interpretation of these data. We used the data collected along with OA gene expression data to expand existing annotations and prioritise the most promising therapeutic targets, which we validated using the latest reported associations. We constructed an online database OATargets to allow researchers to explore the collated data and integrate it with existing OA and skeletal gene expression data. CONCLUSIONS: We present a comprehensive survey and online resource for understanding gene regulation of animal model OA pathogenesis.

Regulation of microRNA-221, -222, -21 and -27 in articular cartilage subjected to abnormal compressive forces.

KEY POINTS: microRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression. miRs are involved in regulating cellular activities in response to mechanical loading in all physiological systems, although it is largely unknown whether this response differs with increasing magnitudes of load. miR-221, miR-222, miR-21-5p and miR-27a-5p were significantly increased in ex vivo cartilage explants subjected to increasing load magnitude and in in vivo joint cartilage exposed to abnormal loading. TIMP3 and CPEB3 are putative miR targets in chondrocytes Identification of mechanically regulated miRs that have potential to impact on tissue homeostasis provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology, in all physiological systems. ABSTRACT: MicroRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression and are involved in mechano-regulation of cellular activities in all physiological systems. It is unknown whether such epigenetic mechanisms are regulated in response to increasing magnitudes of load. The present study investigated mechano-regulation of miRs in articular cartilage subjected to 'physiological' and 'non-physiological' compressive loads in vitro as a model system and validated findings in an in vivo model of abnormal joint loading. Bovine full-depth articular cartilage explants were loaded to 2.5 MPa (physiological) or 7 MPa (non-physiological) (1 Hz, 15 min) and mechanically-regulated miRs identified using next generation sequencing and verified using a quantitative PCR. Downstream targets were verified using miR-specific mimics or inhibitors in conjunction with 3'-UTR luciferase activity assays. A subset of miRs were mechanically-regulated in ex vivo cartilage explants and in vivo joint cartilage. miR-221, miR-222, miR-21-5p and miR-27a-5p were increased and miR-483 levels decreased with increasing load magnitude. Tissue inhibitor of metalloproteinase 3 (TIMP3) and cytoplasmic polyadenylation element binding protein 3 (CPEB3) were identified as putative downstream targets. Our data confirm miR-221 and -222 mechano-regulation and demonstrates novel mechano-regulation of miR-21-5p and miR-27a-5p in ex vivo and in vivo cartilage loading models. TIMP3 and CPEB3 are putative miR targets in chondrocytes. Identification of specific miRs that are regulated by increasing load magnitude, as well as their potential to impact on tissue homeostasis, has direct relevance to other mechano-sensitive physiological systems and provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology.

The role of microRNA-3085 in chondrocyte function.

MicroRNAs have been shown to play a role in cartilage development, homeostasis and breakdown during osteoarthritis. We previously identified miR-3085 in humans as a chondrocyte-selective microRNA, however it could not be detected by Northern blot. The aim of the current study was to prove that miR-3085 is a microRNA and to investigate the function of miR-3085 in signaling pathways relevant to cartilage homeostasis and osteoarthritis. Here, we confirm that miR-3085 is a microRNA and not another class of small RNA using (1) a pre-miR hairpin maturation assay, (2) expression levels in a Dicer null cell line, and (3) Ago2 pulldown. MicroRNA-3085-3p is expressed more highly in micromass than monolayer cultured chondrocytes. Transfection of miR-3085-3p into chondrocytes decreases expression of COL2A1 and ACAN, both of which are validated as direct targets of miR-3085-3p. Interleukin-1 induces the expression of miR-3085-3p, at least in part via NFκB. In a feed-forward mechanism, miR-3085-3p then potentiates NFκB signaling. However, at early time points after transfection, its action appears to be inhibitory. MyD88 has been shown to be a direct target of miR-3085-3p and may be responsible for the early inhibition of NFκB signaling. However, at later time points, MyD88 knockdown remains inhibitory and so other functions of miR-3085-3p are clearly dominant. TGFß1 also induces the expression of miR-3085-3p, but in this instance, it exerts a feedback inhibition on signaling with SMAD3 and SMAD4 shown to be direct targets. This in vitro analysis shows that miR-3085-3p functions in chondrocytes to induce IL-1-signaling, reduce TGFß1 signaling, and inhibit expression of matrix genes. These data suggest that miR-3085-3p has a role in chondrocyte function and could contribute to the process of osteoarthritis.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs.

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

Combined medial patellofemoral and patellotibial reconstruction with soft tissue fixation in recurrent patellar dislocation.

INTRODUCTION: The medial patellofemoral and patellotibial ligaments (MPFL and MPTL) are the main passive restraints to lateral patellar translation. When nonoperative management of patellofemoral dislocations fails, surgical options can be considered to restore patellofemoral stability. Several reconstruction procedures of the MPFL with semitendinosus, gracilis, quadriceps tendon, and synthetic grafts have been described. No clear superiority of one surgical technique over another is evident. MATERIALS AND METHODS: Patients who suffered at least two documented episodes of unilateral patellar dislocation, confirmed radiographically and at clinical examination, underwent combined MPFL and MPTL reconstruction. Patients were regularly followed-up postoperatively at 2, 4, 8, 12, and 24 weeks, and then annually for a minimum of 2.5 years. Clinical and functional evaluations were performed using the modified Cincinnati rating system and the Kujala score, while anthropometry values including thigh volume and cross-sectional area of the thigh were measured before the operation and at the latest follow-up bilaterally. RESULTS: There were 7 males and 27 females with a mean age of 26.5 ± 10.7 years (range, 13-39 years). The mean follow-up was 3.1 years (range, 2.5-4 years). The mean modified Cincinnati score Increased from 51 ± 22 preoperatively to 90 ± 19 (P = .001). The mean Kujala scores increased from 47 ± 17 preoperatively to 82 ± 17 (P = .02), with no significant differences between patients with or without osteochondral lesions (P ≥ .05), and between male and female patients (P ≥ .08). The Insall-Salvati index was 1.1 preoperatively and remained within normal range (P = .05) at the latest follow-up. CONCLUSION: Combined reconstruction of MPFL and MPTL using an ipsilateral autologous gracilis tendon is satisfactory and effective and can be considered as suitable management option to treat recurrent dislocation of the patella. However, randomized studies are needed to compare different techniques. STUDY DESIGN: Case series.

Interplay between genetics and epigenetics in osteoarthritis.

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

Medical simulation utilization among pediatric anesthesiology fellowship programs.

BACKGROUND: Simulation-based education is a mainstay in education of pediatric anesthesiology trainees. Despite the known benefits, there is variability in its use and availability among various pediatric anesthesiology fellowship programs. AIM: The primary aim was to understand the current state of simulation-based education among pediatric anesthesiology fellowship programs and define barriers that impede the development of an effective simulation program. METHODS: This survey-based, observational study of simulation activities within United States-based pediatric anesthesiology fellowship programs was approved by the Institutional Review Boards (IRB) of the authors' institutions. A 35-question survey was developed in an iterative manner by simulation educators (AA, WW, DY) and a statistician familiar with survey-based research (AN) using research electronic data capture (REDCap) for tool development and data collation. Descriptive and thematic analyses were performed on the quantitative and qualitative responses in the survey, respectively, and were stratified with small, medium, and large fellowship programs. RESULTS: Forty-five of 60 (75%) fellowship programs responded to the survey. The presence of a dedicated simulation program director and number of simulation instructors was positively associated with the size of program and years in operation. Dedicated simulation support was variable across programs and was usually present within the larger programs. A positive association also existed for educational activities among all programs mostly based on size of program and years in operation. Protected time was the most commonly cited barrier to having a comprehensive and sustainable simulation program. There was general agreement for establishing a standardized and shared curriculum among fellowship programs. Approximately 70% of simulation programs had no formal simulation instructor training requirement. CONCLUSIONS: Simulation-based curricula are broadly offered by many fellowship programs. Improved collaboration locally, regionally, and nationally may improve educational opportunities for fellowship programs, particularly the small ones. These efforts may begin with the development of a standardized curriculum and formal instructor training programs.

CRELD2 Is a Novel LRP1 Chaperone That Regulates Noncanonical WNT Signaling in Skeletal Development.

Cysteine-rich with epidermal growth factor (EGF)-like domains 2 (CRELD2) is an endoplasmic reticulum (ER)-resident chaperone highly activated under ER stress in conditions such as chondrodysplasias; however, its role in healthy skeletal development is unknown. We show for the first time that cartilage-specific deletion of Creld2 results in disrupted endochondral ossification and short limbed dwarfism, whereas deletion of Creld2 in bone results in osteopenia, with a low bone density and altered trabecular architecture. Our study provides the first evidence that CRELD2 promotes the differentiation and maturation of skeletal cells by modulating noncanonical WNT4 signaling regulated by p38 MAPK. Furthermore, we show that CRELD2 is a novel chaperone for the receptor low-density lipoprotein receptor-related protein 1 (LRP1), promoting its transport to the cell surface, and that LRP1 directly regulates WNT4 expression in chondrocytes through TGF-ß1 signaling. Therefore, our data provide a novel link between an ER-resident chaperone and the essential WNT signaling pathways active during skeletal differentiation that could be applicable in other WNT-responsive tissues. © 2020 American Society for Bone and Mineral Research. © 2020 The Authors. Journal of Bone and Mineral Research published by American Society for Bone and Mineral Research..

Histone ChIP-Seq identifies differential enhancer usage during chondrogenesis as critical for defining cell-type specificity.

Epigenetic mechanisms are known to regulate gene expression during chondrogenesis. In this study, we have characterized the epigenome during the in vitro differentiation of human mesenchymal stem cells (hMSCs) into chondrocytes. Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) was used to assess a range of N-terminal posttranscriptional modifications (marks) to histone H3 lysines (H3K4me3, H3K4me1, H3K27ac, H3K27me3, and H3K36me3) in both hMSCs and differentiated chondrocytes. Chromatin states were characterized using histone ChIP-seq and cis-regulatory elements were identified in chondrocytes. Chondrocyte enhancers were associated with chondrogenesis-related gene ontology (GO) terms. In silico analysis and integration of DNA methylation data with chondrogenesis chromatin states revealed that enhancers marked by histone marks H3K4me1 and H3K27ac were de-methylated during in vitro chondrogenesis. Similarity analysis between hMSC and chondrocyte chromatin states defined in this study with epigenomes of cell-types defined by the Roadmap Epigenomics project revealed that enhancers are more distinct between cell-types compared to other chromatin states. Motif analysis revealed that the transcription factor SOX9 is enriched in chondrocyte enhancers. Luciferase reporter assays confirmed that chondrocyte enhancers characterized in this study exhibited enhancer activity which may be modulated by DNA methylation and SOX9 overexpression. Altogether, these integrated data illustrate the cross-talk between different epigenetic mechanisms during chondrocyte differentiation.

DNA hypomethylation during MSC chondrogenesis occurs predominantly at enhancer regions.

Regulation of transcription occurs in a cell type specific manner orchestrated by epigenetic mechanisms including DNA methylation. Methylation changes may also play a key role in lineage specification during stem cell differentiation. To further our understanding of epigenetic regulation in chondrocytes we characterised the DNA methylation changes during chondrogenesis of mesenchymal stem cells (MSCs) by Infinium 450 K methylation array. Significant DNA hypomethylation was identified during chondrogenic differentiation including changes at many key cartilage gene loci. Integration with chondrogenesis gene expression data revealed an enrichment of significant CpGs in upregulated genes, while characterisation of significant CpG loci indicated their predominant localisation to enhancer regions. Comparison with methylation profiles of other tissues, including healthy and diseased adult cartilage, identified chondrocyte-specific regions of hypomethylation and the overlap with differentially methylated CpGs in osteoarthritis. Taken together we have associated DNA methylation levels with the chondrocyte phenotype. The consequences of which has potential to improve cartilage generation for tissue engineering purposes and also to provide context for observed methylation changes in cartilage diseases such as osteoarthritis.

Correlation of Infinium HumanMethylation450K and MethylationEPIC BeadChip arrays in cartilage.

DNA methylation of CpG sites is commonly measured using Illumina Infinium BeadChip platforms. The Infinium MethylationEPIC array has replaced the Infinium Methylation450K array. The two arrays use the same technology, with the EPIC array assaying almost double the number of sites than the 450K array. In this study, we compare DNA methylation values of shared CpGs of the same human cartilage samples assayed using both platforms. DNA methylation was measured in 21 human cartilage samples using the both 450K and EPIC arrays. Additional matched 450K and EPIC data in whole tumour and whole blood were downloaded from GEO GSE92580 and GSE86833, respectively. Data were processed using the Bioconductor package Minfi. DNA methylation of six CpG sites was validated for the same 21 cartilage samples by pyrosequencing. In cartilage samples, overall sample correlations of methylation values between arrays were high (Pearson's r > 0.96). However, 50.5% of CpG sites showed poor correlation (r < 0.2) between arrays. Sites with limited variance and with either very high or very low methylation levels in cartilage exhibited lower correlation values, corroborating prior studies in whole blood. Bisulphite pyrosequencing did not highlight one array as generating more accurate methylation values. For a specific CpG site, the array methylation correlation coefficient differed between cartilage, tumour, and whole blood, reflecting the difference in methylation variance between cell types. Researchers should be cautious when analysing methylation of CpG sites that show low methylation variance within the cell type of interest, regardless of the method used to assay methylation.